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ABSTRACT 


This thesis presents a Markov model, which, given an indirect fire weapon system's 
parameters, yields measures of the weapon’s effectiveness in providing fire support to a 
maneuver element. These parameters may be determined for a variety of different 
scenarios. Any indirect fire weapon system may be a candidate for evaluation. This 
model may be used in comparing alternative weapon systems for the role of direct 
support of a Marine Corps infantry battalion. The issue of light gun vs. heavy gun was 
the impetus for the study. The thesis also provides insight into the tactic of frequently 
moving an indirect fire weapon to avoid enemy detection, and possible subsequent 
attack. 
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I. INTRODUCTION 


Field Artillery is known as the King of Battle for good reason. Great captains have 
relied upon its awesome firepower to change the course of battle for many years. It 
provides the commander with a force multiplier which is available day and night, has 
no Weather constraints, and fires a wide variety of ammunition, both conventional and 
nuclear. In addition, it does not involve the sending of friendly personnel across enemy 
lines to deliver ordnance. Few other weapon systems can make these claims. 

There is currently much debate in the artillery community as to the type of weapon 
system which is best suited for the direct support mission. Light weight, mobile 
howitzers or mortars provide flexibility and require little airplane or deck space, but do 
they have enough firepower to do the job? Larger, heavier howitzers, such as the M-198, 
have greater lethality, and an increased maximum range, but does this compensate for 
the lesser mobility caused by their size? Quite often it is the weapon system’s 
characteristics which will determine the tactics used, and, ultimately, the level of support 
provided the maneuver element. This thesis will present a model of a firing battery, 
based on the weapon system characteristics of mobility, maximum range, and 
ammunition lethality.! These quantities will be represented by certain parameters in the 


model. The goals of this thesis are: 


1. An analytical model which provides insight as to the effect the choice of weapon 
system has on the level of support provided to the supported maneuver element. 


2. A look at battery survivability as a function of battlefield conditions, artillery 
tactics, and the weapon svstem choice. 


3. A better understanding of the requirements for modeling indirect fire weapon 
systems. 

Chapter 2 will detail the operations of a direct support field artillery battery. 
Chapter 3 will provide a Markovian model for this operation. Chapter 4 provides the 
results of applying the model in a simple example. Chapter 5 investigates the effect the 
decision to move has on the level of support provided by the battery, and on the 
battery's survivability. It will also show that by allowing the convenient assumption of 


exponential sojourn times in the states of the model, very little accuracy is lost, but much 


1 Although it will be presented in terms of artillery, the model is applicable to any indirect fire 
weapon system. 


is gained in mathematical tractability. Chapter 6 gives conclusions and 


recommendations for further study. 


Il. OPERATIONS OF A MARINE CORPS DIRECT SUPPORT 
ARTILLERY BATTERY 


A. THE MISSION 
Marine Corps artillery has the following responsibilities in support of the 
amphibious assault and subsequent operations ashore: 


1. Providing fire in support of maneuver actions and as a part of the overall fire 
support effort to include: 


¢ Close support of maneuver units in combat. 
¢ Counterfire operations against enemy indirect fire systems. 
¢ Deep interdiction fire on enemy command posts, logistical installations, etc. 


2. Provide fire support planning and coordination resources and facilities to all levels 
of force headquarters. [Ref. 1] 


Each artillery unit is also given a tactical mission. It may be: 
e Direct support (DS) 
¢ General support 
e Reinforcing 


¢ General support reinforcing 


Each tactical mission has its own order of priority by which it responds to requests 
for fire support. For example, a battery assigned a reinforcing mission will respond to 
a request for fire from the artillery unit it is reinforcing, before it does so for requests 
from another source. In the case of direct support, the battery will fire the requested fire 
missions from the supported unit before engaging other targets. In addition, for direct 
support, the zone of fire of the DS battery is the zone of action of the supported unit. 
This provides the supported unit commander with a dedicated fire support asset which 
he can base his maneuver plans upon. While any of these tactical artillery missions can 
be critical to the overall mission success, the area where there is the most debate, as to 


the appropriate weapon system to employ, is in the direct support role. 


B. BATTERY ORGANIZATION 

A Marine DS battery is composed of two firing platoons, and a headquarters 
platoon. Each firing platoon has four howitzers. The battery is commanded by a 
captain, and each platoon by a lieutenant. The battery is one of three firing batteries in 
a direct support artillery battalion. The battalion is commanded by a lieutenant colonel. 
The battalion provides direct support to an infantry regiment, while each of the batteries 
support one of the infantry regiment’s battalions. 


C. TACTICS 
The firing battery may be brought into the battle by air, from the sea, or over land. 
The method used is determined by the tactical situation, terrain, transportation 
available, and the speed with which the fire support is needed. Once joined in battle, the 
most common method of moving the battery is over land. In the case of towed artillery, 
the towing is done by a “prime mover” suited for the load the weapons require. 
The battery will occupy a firing position until an event occurs causing it to move. 
These events can be: 
e The battery comes under attack. 
e The current position does not allow the battery to support its maneuver element. 


e The commander feels the threat of attack due to the time spent in the position, or 
the amount of firing done from the position, has reached an unacceptable level. 
When any of the above events occur, the battery will displace to a new firing position, 
provided the force commander does not feel the current need for fire support outweighs 
the advantage the movement brings. 

Determining the point at which the threat of attack is unacceptable is a difficult 
task. In many cases, a battery will have no indication that their position has been 
compromised until it comes under attack. In the past, rules concerning the time to move 
have been subjectively set by the commander. For example, a rule might be “displace 
the battery when you have been in position four hours, or after firing ten missions, which 
ever comes first”. 

The positions the battery occupy are, by doctrine, about one-third the maximum 
range of the weapon system behind the Foward Edge of the Battle Area (FEBA). The 
two firing platoons will occupy separate positions, typically | km. apart. The howitzers 
set up in the position to take advantage of cover and concealment present, and while 


not firing, work on improving it. 


D. TARGETS 

Targets may be acquired by many different means. The most common source of 
targets for a direct support artillery battery is by a Foward Observer (FO). FOs are 
battery personnel attached to the supported maneuver element. Their purpose is to 
detect and request fire on targets in the maneuver element’s area of operation. They also 
help the maneuver element commander plan fire support. Targets may also be detected 
by countermortar/counter battery radar, aerial observers, and intelligence sources 
external to the artillery unit. 

Targets are classified as to their degree of prearrangement. They may either be 
planned, or targets of opportunity. Planned targets have a high degree of 
prearrangement. Their firing data is prepared in advance, and they are either fired 
according to schedule, or on call. Targets of opportunity do not have firing data 
prepared in advance, hence a low degree of prearrangement. Because of this, it takes 
longer to engage a target of opportunity than a preplanned target. Targets are also 


assigned priorities, based on their potential threat to the supported unit. 


E. THE FIRE MISSION 

A fire mission is the act of bringing artillery rounds on a selected target. The 
processing of a fire mission is a complicated evolution, with many tasks to be 
accomplished concurrently. To begin, a FO may spot a target. He will then request fire 
on it by sending a “call for fire” to the battery Fire Direction Center (FDC). In this 
message he will provide such information as: the type of mission, the target location and 
description, and his preferred method of engaging the target. 

The type of mission will usually be adjust fire, fire for effect, immediate 
smoke]suppression, or suppression of enemy air defense. Adjust fire means that one 
round of ammunition is fired by the battery, and a correction to bring the point of 
impact closer to the target is sent by the FO to the FDC. This continues until the last 
impact is within a specific distance from the target, based on the ammunition’s lethality. 
At this point the entire battery engages the target. Fire for effect is used against targets 
whose location is accurately known, or when surprise or speed is a critical factor. In this 
case the entire battery will engage the target on the first volley. Immediate 
smoke]suppression is used if the supported unit is under attack, and speed of response is 


more critical than accuracy. 


Suppression of enemy air defense is used to suppress enemy anti-aircraft weapons while 
friendly aircraft are operating in the area. 

The FDC will simultaneously determine whether the battery is capable of engaging 
the target (i.e., a mission of higher priority is not being processed, the target is within 
range, and the ammunition is available), and if the target is in a restricted fire zone. The 
artillery battalion FDC and the Fire Support Coordination Center will monitor this 
message and determine if the target may be engaged more effectively by another fire 
support asset (e.g., naval gunfire, etc.), or, whether firing at this target will adversely 
effect other maneuver elements. If they are to allow the battery to fire the mission, 
nothing will be said. “Silence is consent” is the rule followed. 

During this time the firing solution for the howitzers is computed. These data 
include which guns will fire, how many rounds they will fire, the type of round and fuze, 
the powder charge, deflection (direction), quadrant elevation, and the method of control. 
This information is transmitted to the gun sections, who then prepare their weapon for 
firing. Firing commences according to the controls established by the FDC. 

When the determined number of rounds have been fired, the FDC will notify the 
FO of “rounds complete”. The FO will either request additional rounds if the effect on 
the target has not been achieved, shift to another target if one is available, or end the 
mission. When the mission is ended, the FO will notify the FDC of his assesment of 
damage to the target. 

The method of attacking a target is influenced by the results desired. This is usually 
to suppress, neutralize, or destroy a target. Suppressive fire is delivered to deny the 
enemy the opportunity to fire his weapon or maneuver freely. Accuracy is not as 
important as is the frequency with which the rounds impact. Neutralization fire 
attempts to degrade the enemy’s combat efficiency. Surprise of the impacting rounds, 
and the ability to mass fire from several batteries, can be critical to the success of 
neutralization fires. Destructive fire attempts to make a target permanently ineffective. 
Accuracy and round lethality are important factors in a target destruction mission. 

Although the fire support operation is complex, it is carried out quite efficiently. 
But, it can be improved. Having the proper weapon system for the job makes it much 
easier, and also more effective. Employing the tactics best suited for the weapon system 
used will enhance the support provided, while increasing the battery’s survivability. 


I. THE MODEL 


A. INTRODUCTION 


The Model is formulated so that it can represent a wide range of weapon types, from 


mortars to heavier self-propelled artillery. This was accomplished through the varying 


of certain model parameters. Different combat conditions, from low-intensity conflict 


through high-intensity conflict, can also be modeled by varying parameters in the 


model. 


Although a direct support artillery battery currently consists of two firing 


platoons, the model will consider it as an entity, to allow more flexibility in the 


comparison of alternative weapons and scenarios. 


B. MODEL STATE SPACE 


The battery will be modeled as an irreducible, recurrent Markov Chain, with state 


space as follows: 


1B 


t 


Q; 


The battery is in a firing position, available to support the maneuver 
element, and has fired i volleys since occupying this position. It has not 
been detected by the enemy. 


The battery is in a firing position, available to support the maneuver 
element, and has fired i vollevs since occupying this position. It has been 
detected bv the enemy, but is unaware of it. 


The battery is in a firing position. It has been detected by the enemy, and 
it is aware of this. Because it is preparing to leave the position, it 1s not 
available to support the maneuver element. 


The battery is displacing to a new firing position. It is not available to 
support the maneuver element. 


The battery is “down” after sustaining an attack. It is not available to 
support the maneuver element. 


C. MODEL PARAMETERS 


The following are the parameters of the model. 


Number of different target types the battery may engage. 
Mission (target) arrival rate of type j targets. j= 1,2,...,J 
Number of volleys to fire at target of type j. ja 2,8 
Prob(detection from firing the i-th single round or volley) 


Rate of response of the enemy once the battery is detected; 1/p is the mean 
time until the enemy attacks the battery, after detection. 


Rate at which the battery repairs itself after sustaining an attack; 
1/6 is the mean down time. 


Rate at which moves, between firing positions, are completed; I/p is the 
mean time of a move. 


Rate at which the battery is attacked while moving; @ dt is the 
(approximate) probability that the battery is attacked at time t after the 
move has begun. 


Rate at which the battery determines its firing has caused detection. 


Rate at which the battery leaves a position, once the decision to move has 
been made. 


Number of rounds/vollevs at which the battery will displace; k is a decision 
variable, the value of which is to be selected. 


Table 1 on page 10 gives an interpretation of these parameters as a function of the 


weapon system and/or the scenario. Appendix A details additional considerations a user 


should make when estimating these parameters. 


D. ASSUMPTIONS 


The model makes the following assumptions: 


I. Each time a battery occupies a new firing position, it is in an undetected state. 


Zale 
countermortar/counterbattery radar. The detection of one round allows the 
accurate locating of the battery. In addition: 


primary source the enemy has for detecting the battery is 


Prob(detecting a single round) = Prob(detecting a volley fired simultaneously) 


Other methods employed by the enemy for locating the battery also rely heavily 
on the batterv’s firing for detection. 


. The probability of detection is not a function of the number of rounds fired 
previously without detection, i.e., 


Pi = Pp i 'e Me, one k. 


. Missions (targets) of different types arrive according to a homogeneous Poisson 
process. The battery will fire a set number of volleys at each target type. 


. The battery will displace to a new firing position after firing a predetermined 
number of rounds/volleys, if it becomes aware of its own detection; see next; or if 
needed to maintain coverage of its supported element’s zone of action. 


. If detected, there is a chance that the battery will become aware of this before it is 
attacked. It will then attempt to leave the position before the attack begins. 


. The battery cannot accept fire missions during its movement between firing 
positions. 


. If attacked in position, the battery will suffer casualties and battle damage, causing 
it to be unavailable to support the maneuver element for a period of time, or 
forever. 


Table 1. PHYSICAL INTERPRETATION OF THE MODEL PARAMETERS 


high values of the parameter low values of the parameter 


Weapon system capable of engaging a Weapon system capable of engagin 
variety of different target types due to its | limited number of different target t 
ammunition and its maximum range. because of ammunition and maxim 
range constraints. 













High intensity conflict. low intensity conflict. 





Large weapon system maximum range. 


_ Be Small, less lethal rounds. Large, more lethal rounds. 
J=1,2.0.. 


An opponent with a very capable An opponent with a poor counterba 


Small weapon system maximum ré 


counterbattery locating system. locating system. 


Very capable enemy. Relatively incapable enemy. 


Battery located near the foward edge Battery located far from the fow 
of the battle area. edge of the battle area. 


Battery one of few targets. Battery one of many targets. 


Enemy response capacity limited. Enemy response means capabl 


Battery in a fortified position. Battery in a vulnerable positio: 


Mobile weapon system. Easily Less mobile weapon system. Diffict 
transported over the terrain. transport. 
Enemy very capable of detecting and Enemy with a lesser capability « 
attacking a moving convoy. detecting and attacking a movin. 
convoy. 
Friendly means of detecting an Friendly means of detecting an 
impending attack is very good. impending attack is poor. 





| a Weapon system able to displace rapidly Weapon system slow to displace 


E. MODEL FORMULATION 
Let: 

P(t) = Prob(being in state P, at time t). 
Q(t) = Prob(being in state Q, at time t). 
R(t) = Prob(being in state R, at time t). 
U(t) = Prob(being in state U at time t). 
D(t) = Prob(being in state D at time t). 
p, = 1-(1—-pys 

q = 1-B, 


The rate equations are: 


J 
PL 1 — > Ade] + Ulude + D(Nsar 


Po(t + at) = 
j=! 
J i-1 J 
P(t+dt) = rio] | z Saat | + DP) x >, Agar i= eee 
j=) i=0 j=) 
sf btn; =i 
k—1 J 
Pred) = PiQli-ad] + PM) x d, Aga 
i=0 j=l 
st den, 2k 
Q(t) = 0 


J 
Oft+d) = Ojt) [1 Diya x= pan xc — rat 4 
J 


=I 


i=] df i-] df 
| >, x > uo | i px x > ua po eal 
i=0 j=l = jal 


5.0 Ltn, =f Sn Ltn, =i 


1] 


Ot+dt) = Of) [(1- adt)x(1—padr)x(1—ydn)] + 


k-1 df k-1 J 
SP) x >. Aydt + DO) x DY Ade 
Sin, 2k So +n; 2k 


R(t+d) = RX) C(1—pd)x(1—ad)] + Old yat pal gel 


k 
Urtd) = Ut) [(1—udt)x(1-@ad)] + Pradt + O,(t)adt + D Rit) a de 
i=] 


D(t + dt) 


k k 
DVL 1—Sdt] + DOW) pdr + UM Od + ) RW) pat 


i=) i=] 


These equations may be solved for P(z), Ot), U(t), R(t), and D(t) [Ref. 2]. Because 
the Markov chain is irreducible and positive recurrent, the limiting distributions exist. 
Solving for them will allow the model to predict the proportion of time the battery will 


spend in a state where it can provide fire support to the maneuver element. 


Lee 
N(z) = lim P{t) 
Q(i) = lim Q(1) 
Qi) = lim Rit) 
U = lim Ut) 
D = lim Dt 


The steady state equations are: 


df J 
by 4. D0 | M0) = »xU+6D 


j=) j=) 


df ff df 
ban Ya) Mi) = » 4,4 MG —n) Poe 
j=1 j=1 j=i 


Gye i-n, 20 
k=-1 
; J 
mA) = +x y y yan 
Pa a 
Sut +n; 2k 
J J J 
= jai =a 
(-n)20 (i-n)21 


k-1 
aj es 
y A,pjT(k—n) + » ys 1, (i) 
Z ry 


j=! 
S.t: i-n, >0 = s.cl+n, 2k 

0) k = tae yi 

© (oe ceca) 
where: 
©(0) = 0 

‘ iv : ; 
Q(z) = oa x O(2) r= ee 


k 
y = —<— x LT) + OW + > 20] 


i=] 


k k 
i a xu + £ x Yo + 4 x Yaw 
l=1 i=1 
For given values of J and n, ,these equations can be solved as a function of I1(0). 
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Let: 
T1(i)’ 


MI (i) + TI (0) 
O(i)’ = O() + TT (0) 
Qi)! = QQ) + TT (0) 
joy 10) 3 1 
U'e=U+T1 (0) 


The sum of the limiting probabilities, over all states, must be one, therefore: 
k k k 

10 + )o@ + YaH + U+DdD=1 

i=0 i=} i=1 


Then: 
k k k 

T1(0) + > Ty x (0) + >. 900! x TNO) + > 20) x (0) + D’x HO) + U xO) = 1. 
i=i i=| i=l 

Solving for I1(0) : 


] 


M10) = k k k 
1+ Ty + Ow + ON + D4 U 
=i = =|) 


Let: m(k) = The long run proportion of time the battery spends in a state in which 


support can be provided to the maneuver element. Then: 


k kK 
mk) = > m9 + Yo 
i= 


i=0 


k k 
=1(0) + )-N1’ x M0) + ) ea x MO) 
imj i=] 


14 


ee Sn a Yew 
f= f=] 


k k k 
1+) ny +) ow+ Yow+ osu 
He) =) =) 


The value of m(k) can be determined for each weapon system, allowing us to 


compare the percentage of total time each would be available, as a function of k. 


IV. AN EXAMPLE 


A. ALGEBRAIC SOLUTION 
This chapter will provide an example of the model solving process. The parameters 
used are arbitrary. To simplify calculations, let: 


j=1 
n=n=1 
4,=2 

2 =P 

Gq, — d= lp 


This is equivalent to there being one type of target available for the battery to 
engage, and the firing of one volley will have the desired effect. A diagram of the model 
is in Figure 1 on page 17. 

The steady state equations are: 


For IT(i): 
J 
dae + dyn] me @= Dd) A4qn-n) F= 12k — I 
Sut a 
Which now simplifies to: 
Axl) = 4g x Ne-)) f= 12,.,.4 =a 
Solving for IT(/) : 
Ni) = ¢ x MG-) i= 1,2,..,4 -1 
This is solved in terms of II(0) : 
ni) = gq’ x M(0) f= 1,2). 


Figure I, 
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Pore: 


k-1 
df 
l y ee 
mk) = +x) YY jg nw 
= 
i=0 


SL itn, >k 


Which now simplifies to: 


Mk) = gL x M(k~1) 


This is solved in terms of I1(0) : 





A q* 
mk) = —— x 1(0) (4.1) 
For O(i) : 
@(0) = 0 
J J J 
j=) j=) j=l 
5.0 ((—n)20 S$. ((-n)21 


Which now simplifies to: 
(y+p+4)O0K%) = ApN-}) + AO00-) i =122=...k a 


This is solved in terms of I1(0) : 


i ! 
i Gs co si 
OQ) = pxCx oan |x T1(0) Pees = | 
Where: 
= | 
om yti+op 


For ©(k): 


k-1 
J Jf 
> = 


j= 

5.0 in, 20 7 5.0 itn, 2k 
og) = 
(ep 7c) 
Which now simplifies to: 


A pM(k—1) + 2@(k—-1) 


ae (are ae 62) 


This is solved in terms of T1(0) : 


; -1 
faor (7): 
oe = 2 
QQ(i) cE - mn x O(2) i 
ForU: 





k 
To o : 
u= 25 x [nw + ow + S20 | 


Where: 





k k K 
ypc | CeCe 
218 ererieeesrey PES =| ary Om 


And I(k) and ©(&) are as determined in equations 4.1 and 4.2, 
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For D: 
a 6 f= 
De 52,00 + 2xu+ 5% 20 


Where: 


: (a+) ¥ 
200 = x) 200 
i=1 t=) 


Since the long run proportion of time spent in all the states must sum to one, 


k-1 


k—1 k 
yng + 1k) + OM + OH + DOH + U4 D=1 
i=] 


i=0 i=1 





Where: 

oun k-1 =. 

dn@ = ya x 110) = — |» T1(0) 
i=0 = gi 

and: 


(i) »4 iG x =e) x (0) 
zl 5 . q—¢ 


i=] 
pc i= Gre c 
“aes «| l—4q - $=€ xn0| 
Therefore, letting: 
T1(i)’ = TI (i) + M1 (0), 
O(i)’ ©) + TI (0), 
Qay’ = Qi) + TH (0), 


D = p= iO), 
U = U=T 0), 
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then: 


k=l k-1 k 
ae = dng + m6 + >" + OK) + DQ" + Dr + U. 
i=0 (=I l=] 


The long run proportion of time the battery is available to support the maneuver 


element, m(k), can be determined for different values of k. 


k—1 k-1 
m(k) = ao x ; DM + Tk) + )O() + eco 


i=) 


B. NUMERIC SOLUTION 


AS an example, let: 


p = 0.05 
q = 0.95 
A= 4. 
p = 2. 
6 = 0.1 
B= 1. 
¢= 0.1 
ao = 6. 


Solutions of m(k) for k from | to 15 are given in table Table 2 on page 22 
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Table 2. MODEL SOLUTION FOR 
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V. DETERMINING THE OPTIMUM MOVE STRATEGY 


The number of rounds or volleys, k, is that which, when fired from the same firing 
position, will cause the battery to move to another position. It may be anticipated that 
if the battery stays in place for too long, measured in the number of rounds fired, it is 
discovered, and attacked, and hence out of action for a prolonged time. If it shoots and 
moves too soon its effectiveness is reduced by being unavailable during many moves. 
Some compromise may be advisable. 

Typically, a firing battery will occupy a position for a predetermined time, or until 
it has fired a predetermined number of rounds. It is assumed that these are the main 
factors in the enemy locating, and subsequently, attacking the battery. Many types of 
locating systems are in use. These include radar, sound and flash ranging, aerial 
observation (visual and IR), as well as human reconnaissance. While not firing, the 
battery has a number of methods to avoid detection, including natural camouflage and 
IR reflecting nets. In most cases it is the battery’s firing which gives its position away. 
It was for this reason that the model was based on the number of rounds fired, for 
determining when to move, and not the time in position. 

In the initial construction of the model, it was thought that the value of k used in 
the comparison of weapon systems would be that k which produced the maximum value 
of m(k), i.e., 


ko = max { m(K) }. 


This implies the commander uses an optimal strategy for moving his firing unit. The 
values of m(k) would be found using the appropiate parameters for the weapon system 
and scenario considered. Graphs of m(k) vs k are in Figure 2 on page 24 ,for arbitrary 
sets of parameters.2 

The most interesting result of this facet of the analysis was that for almost all 


reasonable combinations of parameters a pronounced optimum, not at an extreme point, 


2 In each graph the parameters are as follows: 
1=4,p=0.05,n=1,p =2,0=0.1,6=0.1,y =0.1,@ =6 ; except: 1) The upper left graph, the 
solid line is for p = 0.01, and the dotted line is for p = 0.20; 2) The upper right graph, the solid 
line is for 6 = 0.01 and the dotted line is for 6 = 1 ; 3) The lower left graph, p = 0.20 and, the solid 
line is for 9 = 0.01 and and the dotted line is for 6 = 1 ; 4) The lower right graph, the solid line is 
for y = 0.05 and the dotted line is for y = 2. 
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Figure 2. mj(k) vs. k 
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rarely exists. Quite often, when a non-extreme optimal point exists, as in the dotted line 
in the upper left graph of Figure 2 on page 24, the difference between the m( k’ ) value 
obtained at the optimal point and that at very large values of k is very small. Using this 
method would then result in the decision being either: move after every mission, or do 
not base your decision to move on the number of rounds fired, in almost every case. It 
is of importance to note that the second instance, when m(k) is constantly increasing, 
does not that suggest that one should never leave the firing position. What it suggests 
is not to use the number of rounds fired from that position as the basis for decision. 
Instead, other factors, such as indications of an impending attack, time in position, or 
the need to maintain coverage of the supported unit’s zone of action should be used as 
the basis for the decision. 

When an optimal value exists between extreme points it is not always clear that this 
should be the decision rule. Since the model does not account for additional costs 
associated with moving the battery, such as the human fatigue and equipment wear, the 
small gain which frequent moving sometimes brings may be negated by these costs. 
When the optimum is significantly larger than the values at large k, then the rule would 
be of value, though this rarely occurs. 


These results were viewed very skeptically at first. Artilley has, in the past, always 
used as a major factor in determining when to move a unit the number of rounds fired 
from its position. The determination of this number was usually left to the commander’s 
judgement. It seemed intuitive that firing a great number of rounds from the same 
position would place the battery under increased risk of detection, and subsequent 
attack. In analyzing the model, the memoryless property of the exponentially distributed 
state sojourn times was investigated for its effect on m(k). A Fortran simulation was 
written [Appendix C] in which the Down, Moving, and Response time were made 
nonexponential, with the same mean the Markov model's states had.3 As Figure 3 on 
page 27 shows, the results were very similar, especially at higher values of k. While this 
tended to disprove the possibility of the exponential sojourn times causing the lack of a 
non-extreme maximum it also showed that allowing the convenient assumption of 
exponential sojourn times did not detract from the model’s accuracy very much. This 
assumption made the model much more mathematically tractable. Making the down 


time a constant, equal to the mean of the exponential down time, was also done with the 


3 They were made the sum of five identical exponential distributions. 
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same simulation, with the thought that perhaps it was the variance of the sojourn times 
which caused the given results. The simulation showed that this was not the case. 

Before accepting that the current method of deciding when to move the battery was 
perhaps in need of change, another method of determing k was attempted. By making 
the down state an absorbing state in the Markov model, the expected time the battery 
would be in states where it would provide fire support to the maneuver element before 
absorbtion ( E(T) ) could be determined. An example of this model formulation and 
solving process is detailed in Appendix B. It was felt that possibly this expected time 
would yield a definitive, non-extreme, maximum value, as a function of k. Figure 4 on 
page 28 is a display of two graphs of this expected time vs k. For most realistic 
combinations of parameters, the expected time function behaved as in the left graph of 
Figure 4 on page 28. Rarely was there a maximum not at an extreme point. When it 
did exist, as in the dotted line of the same graph, it was rarely pronounced. The 
difference between the optimal value of the expected time, and the expected time at large 
values of k, was very small.4 

While the results seemed counterintuitive, they were now supported by the two 
models. An explanation was sought as to why this may be the case. One reason for the 
results was that newer countermortar;/counterbattery radars in use today do not need 
several rounds to locate a battery, as they did in the past. If an older system required 
N rounds to zero in on the battery, then moving before firing N rounds from the same 
position would make it impossible for their radar to locate the battery. The new radars 
need only a single round detected for an accurate locating of the battery. Because, like 
a battery firing, they too are subject to being located and attacked, these radars will not 
be active at all times. When they are active, if a battery fires into an area they are 
covering, the battery is detected and located, and presumably targeted for a future 
attack. If they are not active, then the battery has not been located, and it is as if it had 
not fired the last round (for detection purposes). Because it is not known who a future 
advesary will be, or what his radar doctrine will be, the model assumed a constant 
probability, p, of being detected. In actuality, the value of p may be a function of the 
amount of firing a battery is doing. If it is firing frequently, the the radar commander 
will feel that his becoming active will give him a high probability of detecting the battery, 


4 Both graphs have as parameters: 4=4,p=0.05,p =2,u4=1,a=6. The left graph has 
y = 0.05 , and, from lower to higher, @ = 1,0.1,0.01 . The right graph has @ = 0.1 , and, from lower 
to higher, y = 0; & 2: 
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COMPARISON OF MARKOV MODEL AND SIMULATION 
DOTTED LINE IS SIMULATION 


LONG RUN PROB (AVAILABLE TO SUPPORT ) 





4 8 12 
K, THE DECISION VARIABLE 


Figure 3. Comparison of m(k) from the Model and from a Simulation. 
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EXPECTED TIME IN SUPPORT UNTIL ATTACKED VS. K 


EFFECT OF INCREASED RISK DURING MOVEMENT EFFECT OF INCREASED WARNING OF ATTACK 





Figure 4. Expected Time Available to Support Until Attacked vs. k. 
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SO p may increase with the firing rate of the battery. At the same time, though, the 
increased firing may be saturating the capabilities of the radar units, causing a decrease 
in p as the firing rate increases. A constant p reflects the fact that it is not known how 
the parameters will effect the probability of detection. This information should become 
known as a battle progresses. 

Another factor in the shape of the m(k) and E(T) curves is the degree of uncertainty 
that exists with regard to a future attack on the battery. This is shown in the right graph 
of Figure 4 on page 28 . The model allows a probability, ydr , that once detected, the 
battery will know that it is detected in a time, dt. As this value increases, the decision 
to move becomes less a function of the number of rounds fired, and more a function of 
when the battery is informed of an attack. This also depends on the relative response 
rate of the attacking enemy. Knowing of an attack does little good if time isn't available 
to leave the position before it begins. 

To solve the Markov model, and obtain a measure of effectiveness for competing 


Weapon systems, a value of k is required. The calculation 


ki = max { m(k) } 


will provide this value, if it exists, and it is not nearly equal the value of m(k) as k gets 
very large. When k is very small, the result would then be, basically, moving after each 
mission. When m(k) is constantly increasing, another method is needed, because the 
battery will not remain in a position for an indefinite time. In this case, the value of k 


used in the model should be: 


ko = max LEI) }, 
if it exists, Where: 


E(T) = The expected time the battery will be available to support the maneuver 


element, before it is attacked. 


An example of this calculation is in Appendix B. When m(k) and the expected support 


time are constantly increasing, the battery should move: 
* to keep coverage of the supported unit’s zone of action, 


e if attacked, 
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e if warned of a suspected pending attack. 


Then, tet 


M, = Time until battery moves to maintain support. 


M, = Time until battery moves because of an attack. 
M, = Time until battery moves because of a warning of an attack. 
Then: M = Time the battery moves. 


= min{M,,M,,M,} 


Assume that M,,M,, M, are independent and distributed exponentially, with means 
i/u,, 1/u,, 1/u, , respectively. Then: 


F(M) = Prob(M > 1) = eb 7 UstHat Mt} 


It follows that: 


Mea ite th 


EM) = | Ran : 
0 


The expected number of rounds to fire in this time 1s: 


J 
E(k) = E(M)x > Ayn, 
j=) 


which is the value of k to use in comparing weapon systems. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


This chapter presents the conclusions arrived at during the model building effort and 
possible measures of effectiveness. In addition, recommendations for further study and 


a brief summary are given. 


A. CONCLUSIONS 

The initial model, formulated for the comparison of indirect fire weapon systems’ 
effectiveness, was a fairly simple Markov chain. It consisted of k-1 undetected and 
detected states, a moving state, and a down state. The process transitioned from the 
(k-1)-st state to the moving state at the rate at which targets arrived ( 7 ). This 
represented the battery leaving the firing position immediately after firing the k-th round. 
A k-th undetected and detected state was added to the model, with transitions to the 
moving state occuring at a rate different from the target arrival rate. This rate ( « ) was 
representative of the battery’s need to perform certain tasks prior to moving. Since the 
time to complete this preparation is relatively short, the model's long run probabilities 
were changed very little by this modification. Solving for these probabilities was made 
more difficult, so a k-th state was not used in many future models when the goal was to 
determine the effect a secific model change had on the long run probabilites. 

Since there is some risk to the battery in remaining in a firing position too long, but 
there is also a risk associated with moving, it is reasonable to expect that there would 
be an optimal movement decision variable (k). This would yield the largest value for the 
long run proportion of time the battery is available to provide fire support (m(k)). As 
was detailed in Chapter V, this was not always the case. In certain cases, such as when 
the probability of detection (p) was high, or the probability of being attacked while on 
the move was low, the model then determined that to maximize m(k) the battery should 
move after each mission. This conforms to what many in the artillery community feel 
will be the tactic in future conflicts. Artillery will be kept silent until the time when its 
firing will have the greatest impact on the battle. Only then will it fire, and then it will 
move. 

In some instances, a maximum existed for m(k) which was not at an extreme, but 
it was never very pronounced. The difference between the maximum value of m(k), and 
its value at large values of k, was usually small. Since moving the battery entails 


additional costs not measured by this model (e.g., fatigue of the troops and wear on the 
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equipment), the small gain associated with moving frequently would probably not be 
worth these additional costs. 

In some cases, the model predicted that the long run probabilities, m(k), will be 
increasing as k increases. This suggested that the decision to move not be based upon 
the number of rounds fired from the position. Rather, the decision should be based on 
such factors as the need to maintain coverage of the supported unit’s zone of action, or 
because the battery has received some indication that its position has been 
compromised. To determine the effect advanced warning of an attack has, additional k 
states in which the battery has been detected and it knows this, were included in the 
model. These states had little effect on the m(k) values for most combinations of 
parameters. It did have an effect on the expected total time the battery is available to 
provide fire support before being attacked, E(T). As the probability of knowing of an 
attack increased, the maximum in the E(T) curve, if it existed not at an extreme point, 
became less pronounced. This should be expected, as there would be little reason to 
move because of fear of attack, when it is expected that there will be a warning before 
the attack. The presence of a pronounced maximum for E(T) was sensitive to the 
combined effects of some of the parameters. For example, while a high probability of 
knowing of an impending attack tended not to produce a clear maximum, as the average 
response time became shorter, the maximum became more pronounced. This should be 
expected, since knowing of an attack does little good if the response is so quick that the 
battery does not have the time to leave the position before the attack begins. 

The lack of a pronounced maximum in the m(k) and E(T) curves may be due to 
improved systems for locating indirect fire weapons. If so, the tactic of moving based 
upon the number of rounds fired may be in need of change. But, the lack of a 
pronounced maximum may well result from assumptions made in formulating the model. 
Recommendations for further investigation of the effect of certain assumptions on the 
optimal move strategy are presented in the Recommendations section of this chapter. 

1. Measures of Effectiveness from the Model 

Given the value of a set of parameters representing a particular weapon system 
and scenario, and the value for the decision variable, k, determined using methods such 
as those described in Chapter V, the model can be solved for the long run proportion 
of time this weapon system is able to provide fire support to a maneuver element 
(m(k*)). While this in itself is a measure of effectiveness for the support which that 
weapon system provides, other measures may also be determined. A weapon system may 


not be able to engage certain targets because of ammunition or range constraints. While 
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system may be available a greater proportion of time than is a competing system, the 
level of support it provides can be less, if the competing system does not have as great 
an ammunition or range constraint. A possible method of measuring fire support 
effectiveness would be the “gain” a maneuver element commander feels he gets from 
having a particular weapon system supporting him. To measure this gain, a value, V, , 
must be assigned to each of the J target types. This value reflects the advantage a 
maneuver element commander feels is provided him by the battey’s having the desired 


effect on target j. The expected gain per unit of time, could be computed as: 


d 
EL Gain] = m(k’) x )4)x Vy 60 
y=! 


Because J will depend heavily upon the weapon’s maximum range, this will more 
accurately assess the impact of this parameter on the level of support provided. 

This expected gain could also be used in comparing tactics. As an example, a 
commander may locate his firing unit further behind friendly lines, thus increasing its 
security, but decreasing the targets it can engage. By changing the model’s parameters 
to reflect this, a comparison of the expected gain in each case would provide a measure 
of the effectiveness of this tactic. 

The expected total time the battery was available to support the maneuver 
element, until it is attacked, can be determined as is demonstrated in Appendix B. This 
may prove to be a more important measure in certain circumstances, such as when 
replacements for attacked personnel and equipment are not available, or will take an 
extremely long time to acquire. 


Expected costs can also be determined from the model. 


£ 
E({ Cost] = m(k’) x yay ny, 61 
j=l 


is the expected cost, in terms of the number of rounds of ammunition per unit of time 
which is needed to obtain the expected gain. This cost may then be transformed into 


logistic weight, and/or volume, whichever is more appropiate for the analysis. 
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B. RECOMMENDATIONS 

Time did not permit the analysis of several factors impact on the model. These are 
presented as recommendations for further study. 

1. The Movement Decision Strategy. 

The problem of determining the optimal move strategy has been detailed in 
Chapters V and VI. A possible cause for the lack of a more pronounced optimum in the 
m(k) curve may be due to the assumption of a constant probability of detection, p. One 
method to investigate the effect of p on the values obtained for m(k), would be to make 
p an increasing function of the number of rounds fired from the position. This would 
represent an increasing hazard for the battery’s staying in position. Another method 
would be to increase p during periods of relatively high firing rates and decrease p during 
periods of low firing rates. This would model the enemy’s greater willingness to risk 
detection, from radiating his radar, because of his perceived increased probability of 
detecting a fired round. 

If it appears that m(k)’s lack of a pronounced maximum is not due to the model, 
but accurately reflects current counterbattery locating system’s capabilities, then moving 
a battery based on the number of rounds fired from a position may not be a correct 
tactic in many instances. The effect of time in the position, or abilitv to cover the 
supported maneuver element’s zone of action, or even the level of intelligence the unit 
can expect to receive, may be the driving force in determining when to move. 

2. Model Refinement. 

The factors used in creating this model are not exhaustive. They were chosen 
because they were felt to be those with the greatest impact on the effectiveness of fire 
support a weapon system could provide. By modifving existing parameters, and 
including other factors, the model will more closely resemble the operations of a firing 
battery. 

A possible modification would be to let target arrivals be a function of time, 
such as in a non-homogeneous poisson process. Perhaps this function of target arrivals 
may be represented by a Lanchester-type combat model. 

The model can be ammended to reflect the effect of mechanical reliability on 
availability. This could be done as a function of the time spent moving, and the number 
of rounds fired, to model the effect these events have on equipment availability. 

When target arrivals are very fast, or the weapon system’s rate of fire is slow, 
or the number of rounds/volleys required to effect the target is large, there are bound to 


be targets arriving before previous missions have been completed. Since the rate of fire 
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and number of rounds required for effect are specified for a specific weapon system and 
target, the model can be changed to allow targets to join a queue. The order of firing, 
of the targets in the queue, will depend on the the priority assigned to each target. This 
will give an additional measure of a weapon system’s parameters effect on fire support. 
(e.g., the average time a target spends in a queue, or the average length of a queue). 
3. Other Uses for the Model 

This model may be of use for other military systems. The radar systems, which 
attempt to locate a target, operate in a fashion similar to that of a firing battery. It 
actively emits radiation for a period, which while it provides benefits (target locating), 
it also makes it vulnerable to being located, and subsequently attacked. Radars usually 
do emit, for a period of time, then relocate. This model may be of help in evaluating 
radar systems and tactics. 


C. SUMMARY 

The purpose of this thesis is to provide those who must choose the artillery weapon 
system which the Corps will fight with in the next conflict, with a tool to better help 
them make this decision. It is not intended to provide a definitive answer to the question 
of light gun vs. heavy gun. It is intended to demonstrate, within the limitations of any 
probabilistic model, the effects of important weapon system parameters and 
environmental conditions on some measures of the level of support a weapon can 
provide. As with any model, the results are only as good as the assumptions used to 
formulate the model and the accuracy of the estimated parameters. Implementing the 
above recommendations should improve the model. While the decision maker must 
estimate the parameters, this also must be done if a simulation were to be used. When 
a decision maker states that he “feels” that a particular weapon system is more suited for 
a task, he is implicitly estimating these parameters also. The model is able to take these 
estimated parameters and show how their combinations effect fire support effectiveness. 
In addition, the effect of changes in parameters are easy to evaluate, especially when 
compared to performing a large scale simulation. 

Deciding which weapon system to arm our forces with must include some attempt 
at predicting the future battlefield. While, predicting the future is, in general, risky it is 
much more so for combat operations. The cost of being incorrect is much higher, and 
the uncertainty greater. Nonetheless, it must be done. This model may be of some help. 
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APPENDIX A. CONSIDERATIONS FOR PARAMETER VALUES. 


This appendix will give a little more insight into the physical considerations which 


should be made, in determining the parameters to be used, when implementing this 


model. 


Parameter 


J 


This parameter represents the number of different targets and missions 
which the firing battery may engage. As such, it reflects the combination 
of a specific target, e.g., a tank platoon, and a specific mission, e.g., “fire for 
effect”. This is necessary because the result of artillery fire will have a 
different “value” to the Commander for different targets. Also, the same 
target, but engaged in a different mission will require the firing of a different 
amount of ammunition. This will lead to a different probability of 
detection, and a different expected cost, in terms of ammunition fired. 


A weapon system’s maximum range greatly effects J. A weapon with a 
relatively long maximum range will be able to engage targets which are out 
of range of a weapon with a shorter maximum range. The variety of 
ammunition a weapon system is capable of firing will also effect J. If there 
are targets which may be ranged by the battery, but whose ammunition has 
very little, if any, effect on the target, then that target should not be counted 
in J. It would be very unlikely that it would be assigned as a target to the 
battery. 


The level of intensity of combat which the maneuver element is engaged in 
will effect J. It is assumed that combat with Warsaw Pact forces will 
provide a wider variety of target types, than would combat with a lesser 
developed nation. An additional concern is the means of locating a target. 
If the scenario depicted has little friendly targeting assets, than having the 
range and ammunition to attack a target is of little value if the target is 
never located. 


Type j targets (j= 1,...,J), arrive at the battery at the rate A,. As such, 
fy 


Z ; 
oy 
v= 


is the probability that the next target to arrive will be of type j. 





This rate is greatly effected by the intensity of the combat. It is also effected 
by the amount of other friendly fire support assets available. The greater 
the number of options available for fire support, the less frequent will be the 
assignment of targets to the battery. 
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Pi 


Target j will have n, rounds/volleys fired at it, and, it is assumed that this 
will have the desired effect on the target. As an example, if current tables 
indicate that 24 rounds of a given caliber must be fired at a target, then 7, 
for an 8 gun battery would equal 3, (if the mission were fire for effect). If 
it Were an adjust fire mission, an additional number of “adjusting rounds” 
should be added to m, The major factors influencing x, are the ammunition’s 
lethality, the weapon’s firing errors, and the target locating errors. 


The parameter p, is the probability that the battery will be located when 
firing the i-th round/volley from the same position. While the model 
assumed this was constant, this need not be the case. The value assigned 
p, Should reflect the firing signature of the weapon, the rate of fire required 
by the combat conditions, and the quantity and quality of the enemy’s 
indirect fire weapon locating systems. A relatively long maximum range can 
contribute to reducing the probability of being detected by non-radar 
methods, if the battery is positioned further behind the Foward Edge of the 
Battle Area (FEBA). 


The response rate of the enemy is modeled as p . As such, 1/ p is the 
average time it takes for the enemy to attack the battery, once it has 
detected it. 


This parameter has several factors affecting it. One is the importance the 
enemy places on attacking the battery. In a high intensity combat scenario, 
a nuclear-capable artillery battery will have a high priority for being 
attacked, once located. In a situation where the use of nuclear weapons is 
extremely unlikely, the fact that the weapon is capable of firing these 
munitions would have little effect. Also, if the battery's fire support is doing 
much harm to the enemy, then neutralizing the battery will assume a greater 
importance for the enemy. 


The greater the number of similar priority targets (such as additional 
batteries) for the enemy to engage, the longer will be the average time until 
a specific battery is attacked. The response means available to the enemy 
is verv important. If it is limited, then the response will probably take 
longer to occur. Air superiority should also be considered, since it will effect 
the number of options available to the enemy in attacking the battery. 


The weapon’s maximum range, and the tactic used to exploit it, contribute 
to p. A long maximum range allows the weapon system to locate further 
behind the FEBA, thus limiting the enemy’s response options, while still 
maintaining coverage of the supported units zone of action. 


The parameter « represents the rate at which the battery leaves a position, 
once it has been determined that it should leave. Thus, the expected time 
required to prepare to move is I/ a. In addition, 


— 
a sp 


a7 


is the probability that the battery will be attacked before it leaves a position, 
if it has been detected and it is aware of this detection. This parameter is 
also effected by the weapons mobility, and the terrain over which it must 
operate. 


The parameter yp is the rate at which moves are completed. The averagetime 
to complete a move (and be prepared to fire at targets) is, then, 1/ 4. The 
larger, more difficult to transport weapons will have a smaller value for p , 
when compared to a lighter, more mobile system. Additional considerations 
are the terrain expected to be traversed, and the effect of weather on the 
terrain. Another concern may be the weapons size if it limits the number 
of firing positions which the battery may occupy. This could cause longer 
movement distances, thus, longer average move times. 


The parameter 6 represents the rate at which the enemy detects and attacks 
the battery when the battery is moving. As such, the probability that the 
battery is attacked before completing the move is: 


aes 
6+ yn 


This parameter will depend on the signature the battery creates while 
moving, as well as the factors of enemy response, such as air superiority and 
available response means. 


The rate at which the battery determines it has been detected is model as the 
parameter y. The model is formulated such that 


Niet 
Pode op 


is the probability that, given the battery has been detected, it learns of this 
prior to it being attacked. The value of y is largely a function of friendly 
capability to gather and pass to the battery, intelligence concerning enemy 
actions. 


The parameter 6 represents the rate at which the battery recovers from an 
attack. The mean time, from attack until the battery is prepared to accept 
fire missions, is 1/6. 


The value of this parameter will depend upon the severity of the attack, and 
the “hardness” of the battery at the time of attack. Of importance also, is 
the ability of the friendly unit to replace presonnel and equipment lost 
because of the attack. 
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APPENDIX B. EXAMPLE OF THE E(T) CALCULATION. 


The calculation of the expected time a firing battery will be available to provide fire 
support until it is attacked is given for a simple example. This example will use k-1 
states instead of k, for the reason given in Chapter VI, and will use the same 
assumptions as the example in Chapter IV, ie., let: 


a 
n=n=1 
A= 
Pomp 
74 


This is equivalent to there being one type of target available for the battery to 


engage, and the firing of one volley will have the desired effect. 


A. MODEL STATE SPACE 


P, The battery is in a firing position, available to support the maneuver 
element, ard has fired i volleys since occupying this position. It has not 
been detected by the enemy. 


Q; The battery is in a firing position, available to support the maneuver 
element, and has fired i volleys since occupving this position. It has been 
detected by the enemy, but is unaware of it. 


R; The battery is in a firing position. It has been detected by the enemy, and 
it is aware of this. Because it is preparing to leave the position, it is not 
available to support the maneuver element. 


U The battery is displacing to a new firing position. It is not available to 
support the maneuver element. 


A The battery is “down” after sustaining an attack. It is not available to 
support the maneuver element. This is an absorbing state. 
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B. MODEL PARAMETERS 


J Number of different target types the battery may engage. 

A Mission (target) arrival rate. 

n Number of volleys to fire at the target. 

p Prob(detection from firing a single round or volley) 

p Rate of response of the enemy once the battery is detected; 1/p is the mean 
time until the enemy attacks the battery, after detection. 

B Rate at which moves, between firing positions, are completed; 1/u is the 
mean time of a move. 

7] Rate at which the battery is attacked while moving; @ dt is the 
(approximate) probability that the battery is attacked at time t after the 
move has begun. 

y Rate at which the battery determines its firing has caused detection. 

o Rate at which the battery leaves a position, once the decision to move has 
been made. 

k Number of rounds/volleys at which the battery will displace; k is a decision 


variable, the value of which is to be selected. 


C. MODEL FORMULATION 
Let: 
P(t) = Prob(being in state P, at time t). 
Q(t) = Prob(being in state Q, at time t). 
R(t) = Prob(being in state R, at time t). 
U(t) = Prob(being in state U at time t). 
A(t) = Prob(being in state A at time t). 


The rate equations are: 


Pxtt+dt) = P(j(1— Adt) + P_\oAqat i= 1,2,....k—1 


Q(4) = o 
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Ofet+dt) = Of) (C1 — Adt\(l1—pad)(l-ydn] + 


P_\()Apat + O,,() ids a 
R(tt+dt) = Ra) C(1—pa)(l—ad)] + O(nyar = eee 
k-1 


Ue+d) = UNE -nat)(1-6a)] + Py (Qddt + Q%(QAd + DY RWad , 


i=] 


k-] k-1 
A(t+d) = A) + )OWpadr + Ried + UlNed. 
i=} =i 


These yield the following first order differential equations: 


P(t) = -AP{t) + P4q $= 1,2,..k—1 

C= s(t y pO) + OL) + Paw@ip $= 12k —1 

R(t) = —(a + p)R) + Ofdy $= 1,2). 
k-1 

Ui = -(4 + 0)U® + Pi + O44 + DR) a 


k-1 k-1 
A) = 0U0) + YOMp + DRO p 
i=] 


i=! 
These differential equations can be solved by applying a Laplace transform of the form: 


x{s) = | coal: 
0 
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This yields the following equations: 


S p(s) a A pis) as P;-\(s) A q i= re: = 
sqf{s) = —(A+yt pals) + gsls)4 + pls)Ap Nereis | 
yy é 
1) > seem at) i122 ko 

J k-1 
— nt & | 
CO Meirerarre 7) (UO) ech(O) gata «Doh 
P k-1 z k-\ P 
as) = & x us) + Dax 2 + Yn x 4. 


i=1 t= 


These equations are now solved in terms of p,(s): 





Aq i : 
PAs) = (= ae 7] Pols) i=1,2,....4— 1 (4.1) 
P| 4q \ fi-c ” 
qs) = pcx( ae =} | a [pots f= 1,2,..4—1 (A) 
where: 
ew Ats 


GiCs etme | 0s) ee 





pyc Aq _\; Lae 
r(s) = ae eee ) 4 Se fhe f=) 2 ike 
1 k-] 
SO Marietta Gels)) = reper (4.3) 


=) 


The values of p,_,(s) and q,_,(s) can be determined from equations A.1 and A.2, and: 
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k-1 


k-} 
v 
> (9) = (eas. en) « ) as) 


i=| 


k-1 


1~( Aq in 
vic “ mig eo) dq (GO ae) 
~ 


(qB—(s+2)) ] (s+) 


i=] 


and: 


es yp, 


eal (=) 
ase £ u(s) + > ads) - as dhs) = 
i=) i=] 


Where u(s) is determined in equation A.3, and 
k=l KI 
Via) + and) rfs) 2 
i=] 


i=1 


are determined in equations A.5, and A.4. 
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Let: 


1 ifthe battery can provide support att 
Zt) = 
0 Otherwise. 


E(T) = Expected total time the battery is available to provide fire support until 


attacked. Then: 


ge) [Pron Zr) = at 


Since: 
k—1 


Prob Z(t) = 1) = 1— UW) — DRM — ACO, 


i=] 


k-1 
EC i (1 - ui) — ) 2 - ay )a 
i=) 


k-1 
lim 4 - u9)— ) {9 - aa 


s30 : 
1 


k-1 k-1 
= lim | debs + af 
i=) 


i=0 


The summations of p{s) and g{s) have been solved in terms of p,(s) . Since the sum of 
the probability of being in any state, over all states must be 1, p,(s) may be determined. 


er 
PAS) = 7S) as) 
q{s)’ = gs) + pols) 


n(s)’ = n(s) = pds) 

u(s)' = u(s) + p(s) 

a(s)’ = a(s) + p(s) 
Then: 


1 


ae a a 
a Yrs) + sx) afs)’ + sx) nfs) + su(s)’ + sa(s)’ 
i=0 i=l i=) 
and, 


k=l k-1 
E(T) him ( ) pis)’ + as’ x pols) \ 
i=l 


i=0 
Solving for the expected total time the battery is available to support the maneuver 
element before absorbtion (attack), yields: 


k-1 k-1 


Yplsy + Yah 
i=0 i=] 


aS LS ee 
sx ) pfs) + sx ) afs) ae sx > 7(s) ssl S) ee Sia( 5) 
i=0 i=) f=1 
ee pA . eae ee 
P gqA+t+yt+p)—A q P l-c 
i i .=—  - — =< ee 
where: 


€ Ve=l 
Oi gk! pe l= (a) 
p= (Be )x(1+ 7 fee a 
py ady pa ie ae = [| 
(0+ eo + erty sort «|e P | ° aa 


and: 
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APPENDIX C. VS FORTRAN SIMULATION CODE 


TART TTR IK TIT LICE FORKED BC PEN BENET FERED IERIE ACRE RITTER RAI REE 


* A PROGRAM TO SIMULATE THE EFFECT OF THE DECISION VARIABLE K 

* (THE NUMBER OF ROUNDS FIRED FROM A POSITION BEFORE MOVING), ON 
* THE LONG RUN PROPORTION OF TIME SPENT IN STATES WHERE FIRE 
*SUPPORT MAY BE PROVIDED. 


* 
* KEY VARIABLES: , 

* K = NUMBER OF ROUNDS FIRED BEFORE MOVING. DECISION VARIABLE 

* P = PROBABILITY OF DETECTION FROM FIRING A SINGLE ROUND/VOLLEY 

* N = NUMBER OF ROUNDS TO FIRE AT A TARGET 

* TCSMO = TIME TO leave position (CSMO). TCSMOU = MEAN TIME TO CSMO 
* TMOVE = TIME TO MOVE. TMOVEU = MEAN TIME TO MOVE. DIST AS SUM OF 
a L EXPONENTIAL IID RANDOM VARIABLES. 

* TTARG = TARGET INTERARRIVAL TIME. TTARGU = MEAN TIME. 

* TRESP = TIME TO RESPOND, FROM DETECTION. TRESPU = MEAN TIME 

* DIST. AS SUM OF L IID EXPONENTIAL RANDOM VARIABLES 

* TDOWN = TIME DOWN. TDOWNU = MEAN TIME DOWN 

oe DIST AS SUM OF L IID EXPONENTIAL RANDOM VARIABLES 

* TIME = TOTAL TIME FOR THE EVENTS 

* BADT = TIME SPENT IN THE MOVING AND DOWN STATES 

* IT = NUMBER OF ITERATIONS 

* 


dedetedetetetetedete fetete teste tedese te tee Fe dete Fe dete hese Fe KFC VT RII ILI TTT TTT TKR 


PARAMETER (K=1, IA=12245, IB=33456, IC=35567, ID=45578, 
1IE=567893, IG=93765, IT=100, L=5) 


REAL TCSMO, TMOVE, TTARG(K), TRESP, TDOWN, T1(L), T2(L), T3(L) 
1,UZERO1(K), TIME(IT), BADT(IT) 


INTEGER DET 
CALL EXCMS('FILEDEF 11 DISK ARTYSIM OUT Al') 


* SET THE PARAMETER VALUES 
P=0. 05 
TCSMOU=10 
TMOVEU=60 
TDOWNU 
TRESPU 
TTARGU 
N=1 
* TMAX IS THE MINIMUM LENGTH OF TIME AN ITERATION WILL RUN 
TMAX = 100000. 


600. 
30. 
15. 


WRITE(11,*) 'K',K 

WRITE(11,*) 'CSMO MEAN' , TCSMOU, "MOVE MEAN' , THOVEU 

WRITE(11,*) 'P',P, ‘DOWN *MEAN' ,TDOWNU, 'RESPONCE MEAN' , TRESPU 
WRITE(11,*) 'TARGET ARRIV MEAN' ,TTARGU,'MAX TIME' \TMAX, 'IT' , IT 
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DO 888, J=1, IT 


TIME(J) 
BADT( J) 


0 
0 


* GENERATE THE RANDOM NUMBERS 


5 


CALL LEXPN(IA, TCSMO, 1, 1, 0) 
CALL LEXPN( 1B, 11) Ueno) 
CALL LEXPN(IC, TTARG, K, 1, 0) 
CALL IDXPNCID, [2,0 eeu) 
CALL EBXPNCIE, T3, 0, sro) 
CALL LRND(IG, UZERO1, K, 1, 0) 
TMOVE=0 

TRESP=0 

TDOWN=0 

DO 3, I=1,L 


TMOVE=TMOVE + T1(1) 
TRESP=TRESP + T2(1) 
TDOWN=TDOWN + T3(I) 


DO 6, f= ik 
TTARG(I) = TTARGU*TTARG( 1) 


‘TCSMO = TCSMOU*TCSMO 

TMOVE = TMOVEU*TMOVE/L 

TRESP = TRESPU*TRESP/L 

TDOWN = TDOWNU*TDOWN/L 
DET = 0 


PD = 1-(1-P)**N 


%* DETERMINE IF AND WHEN FIRING CAUSES A DETECTION 


10 


DO 10,51 = K/N, 11 
IF(UZERO1(1). LE. PD) THEN 


DET = I 


END IF 
CONTINUE 


* DET = 0 MEANS NO DETECTION OCCURRED FROM THIS FIRING POSITION 
* THE FIRING UNIT WILL THEN DISPLACE. 
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IF(DET. EQ. 0) THEN 


DO 20, I= 1, K/N 
TIME(J) = TIME(J) + TTARG(I) 


TIME(J) = TIME(J) + TCSMO + TMOVE 
BADT(J) = BADT(J) + TMOVE 
GO TO 555 

END IF 
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* DETECTION OCCURS DURING FIRING 
TECDET.NE- ©) THEN 


RTIME = 0 
DO 30, I=DET+1, K/N 
30 ROUMET— RIIME + TTARGC 1) 


RTIME = RTIME + TCSMO 


* DETERMINE IF THE RESPONCE TO THE FIRING OCCURS BEFORE THE FIRING 
* UNIT DISPLACES. 


IF(RTIME. LT. TRESP) THEN 
DO 35, I = 1, DET 
35 TIME(J) = TIME(J) + TTARG(1) 
TIME(J) = TIME(J) + TMOVE + RTIME 
BADT(J) = BADT(J) + TMOVE 
GO TO 555 
END IF 
* RESPONCE HAS OCCURRED BEFORE DISPLACING 
DO 40, I=1, DET 
40 TIME(J) = TIME(J) + TTARG(I) 


TIME(J) = TIME(J) + TRESP + TDOWN 
“BADT(J) = BADT(J) + TDOWN 


END IF 
555 IFC(TIME(J). LT. TMAX) THEN 
GO TO 5 
END IF 
888 CONTINUE 


* STATISTICAL SUMMARY AND OUTPUT 
PGOOD = 0 


DO 666, I=1, IT 
666 PGOOD = PGOOD + (TIME(1)-BADT(I))/TIME(I) 


SUM = 0 


DOmi7e © = 1, IT 
777. SUM = SUM + (((TIME(I)-BADT(I))/TIME(I)) - PGOOD/IT)**2 

VAR = (1/(1T-1. ))*SUM 

SD = VAR**(.5) 


WRITE(11,*) ‘ESTIMATE OF PROPORTION IN SUP. STATE IS',PGOOD/IT 
PRINT *, ‘ESTIMATE OF PROPORTION IN SUP. STATE IS',PGOOD/IT 
WRITE(11,*) ‘ESTIMATE OF VAR OF PROPORTION IN SUP. STATE IS',VAR 
PRINT *, ‘ESTIMATE OF VAR OF PROPORTION IN SUP. STATE IS',VAR 
WRITE(11,*) ‘ESTIMATE OF SD OF PROPORTION IN SUP. STATE IS',SD 
PRINT *, ‘ESTIMATE OF SD OF PROPORTION IN SUP. STATE IS',SD 
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WRITE(11,*) 
END 


! 
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